Suppression of correlated electron escape in double ionization in strong laser fields 
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The effect of the Pauli exclusion principle on double ionization of He atoms by strong, linearly 
polarized laser pulses is analyzed. We show that correlated electron escape, with electron momenta 
symmetric with respect to the field polarization axis, is suppressed if atoms are initially prepared 
in the metastable state 3 S. The effect is a consequence of selection rules for the transition to the 
appropriate outgoing two-electron states. We illustrate the suppression in numerical calculations of 
electron and ion momentum distributions within a reduced dimensionality model. 
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The interaction of atoms with strong laser pulses of 
intensities in the range of 10 14 W/cm 2 may result in the 
production of doubly charged ions. Early experiments 
indicated that the rates for this process are much higher 
than expected on the basis of a single active electron 
model (see [l[ and references therein). The commonly 
accepted explanation is based on the rescattering pro- 
cess, where a temporarily ionized electron is driven back 
to the atom by the electric field p], Q • The acceleration 
by the field provides additional energy which can then 
be transferred to the second electron during the colli- 
sion with the residual ion. Measurements of the doubly 
charged ion and electron momentum distributions, in ad- 
dition to the total ionization yield, [3, HI 0] revealed that 
substantial number of the electrons escapes with equal 
momenta along the field polarization axis. The origin 
of this correlation was analysed as a consequence of the 
Coulomb repulsion between the outgoing electrons, their 
Coulomb attraction to the nucleus and the pulling force 
from the electric field Q. Except for the presence of 
the field during the decay of the complex formed in the 
rescattering this analysis has many parallels to Wanniers 
0, 0] reasoning for double ionization after collision or 
photo excitation. Since in the Wannier case it can be 
shown that the observed cross sections are constrained 
by the symmetries of the problem [H, [^, [l(J [H| , we here 
search for the corresponding effects in strong field double 
ionization. 

Consequences of the exchange symmetry on double 
ionization yield have previously been studied [H, 03 
within the aligned electron model 14j. These authors 



found an reduction in ionization rate in the antisymmet- 
ric state, but overestimation of the Coulomb electron re- 
pulsion in the aligned electron model did not allow them 
to treat the final momentum distributions. These stud- 
ies were motivated by experimental observation of the 
difference in single and double ionization of O2 and N2 
molecules where the initial states are of triplet and sin glet 
character, respectively [15j . Recent experiments [la. Il7|| 
observed differences between O2 and N2 molecules when 
the electron momentum correlations for components par- 



allel to the light polarization axis were measured. 

The constraints on the double ionization cross sec- 
tion arise from the symmetries of the wave functions 
for the outgoing electrons and can be derived from the 
study of the selection rules for transitions to two elec- 
tron continuum states by Maulbetsch and Briggs [TT[ | . 
They analyze transition matrix elements from an ini- 
tial state ip to continuum final states of well defined 
electron momenta ki and k 2 , i.e., (kxk2\T\if>) where T 
is the appropriate transition operator. As r\, r 2 and 
r*i2 = \t\ — r%\ tend to infinity, the continuum states for 
a two electron case of He approach the asymptotic be- 
haviour (fir^l^i^) — * (27r) -3 exp(ifci • r*i + ik2 ■ r 2 +«x), 
where the phase x = 21n(fciri + k\ ■ r\)/ki + 21n(fc2^2 + 
fc2-f 2 )/fc2-ln(fcri2 + fc-ri2)/(2fc) with k = (ki-k 2 )/2 ac- 
counts for the long range of Coulomb potentials. These 
states and the full transition matrix elements can now 
be projected onto a complete set of two electron states 
characterized by total angular momentum L, its projec- 
tion M on a chosen quantization axis z, total spin S 
and any other necessary quantum numbers, collectively 
labelled a, 

(hk 2 \T\ip)= (kifa\aLMS)(aLMS\T\i>). (f) 

uLMS 

Selection rules are now obtained for sets of k\, k%, L, 
M and S for which the overlaps (kifal&LMS) vanish 
identically. Various situations are discussed in [ll|. The 
case most appropriate for correlated double escape with 
equal energy and equal momenta parallel to the field axis 
is their case (I), where fci = k 2 , both momenta lie in a 
plane, i.e., \4>i — <f> 2 \ = ft, and their projections on the 
quantization axis coincide, 0\ = 2 '- then states with 
(S + M) odd do not contribute to the transition matrix 
elements. Therefore, if the initial wave function belongs 
to this subspace and if the evolution operator preserves 
this subspace, these final states will be suppressed. 

For the case of strong field double ionization the tran- 
sition operator T is the evolution operator corresponding 
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to the Hamiltonian (in atomic units) 
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+ F(t)( Zl +z 2 ), (2) 



where F(t) contains the time dependence of a linearly 
polarized laser pulse. It is symmetric under exchange 
of electrons. Therefore, if the initial wave function be- 
longs to the 3 S state of He, the sum in ^ remains re- 
stricted to M = and S = 1, since both the z component 
of the total angular momentum and the spin operators 
commute with the Hamiltonian. As a consequence, the 
matrix element of the transition operator between the ini- 
tial metastable state and continuum states with k\ and 
ki symmetric with respect to the field polarization axis 
vanishes by the above selection rule. Thus, the differen- 
tial cross section for the double ionization process is ex- 
actly zero whenever the momenta of outgoing electrons 
are symmetric with respect to the field axis. Note that in 
experiments starting with the He ground state, i.e. S, 
and in the double ionization of noble gasses the symmet- 
ric escape is not forbidden — correlated electron escapes 
are clearly visible in electron momentum distributions in 
the double ionization of Ne and Ar atoms [J d, 0| • 

On the level of the classical dynamics of the two- 
electron system, the selection rule eliminates a subset 
of phase space that includes the Ci v subspace of sym- 
metric electron motion, which has been shown to cap- 
ture much of the observed momentum distributions [5(. 
In cylindrical coordinates, the C 2v subspace is defined 
by R = Pi = p2, PR = 2p Pl = 2p P2 , Z = zi = z 2 , 
Pz = 2p Zl = 2p Z2 , |0i - 02 1 = I" and = p^ 2 = 0. The 
Hamiltonian ^ in the C 2v subspace reduces to 



h _ Pr+P 2 z 4 

4 Vi? 2 + Z 2 2R 



2ZF{t). (3) 



The potential energy in ([3]) [and obviously also the po- 
tential energy in possesses a saddle located at R = 
R s sm0 s , Z = R s cos8 s , where R s = / y/\F(t)\ and 
8 S = 7r/6 or 5w/6 depending on the sign of F(t). Elec- 
trons moving in the subspace may simultaneously escape 
from the atom by going over the saddle. Trajectories liv- 
ing in the full space that approach the saddle sufficiently 
symmetrically may lead to simultaneous escape of elec- 
trons with highly correlated final momenta [5[. For the 
static electric field (i.e. F(t) =const) the stability analy- 
sis of the saddle allows one to estimate the energy depen- 
dence of the classical cross section for correlated double 
escape for energy close to the saddle energy, i.e. to obtain 
a counterpart of the Wannier threshold law in the pres- 
ence of the electric field [Hj]. Note that for the energy 
equal to the threshold value the only trajectory leading 
to the double escape corresponds to two electrons moving 
symmetrically in the same direction along the field axis. 
This is very different from the original Wannier problem 
where the electrons escape symmetrically but in oppo- 
site directions. The analysis of the quantum selection 



rule for double ionization of the metastable He atoms by 
linearly polarized laser pulses implies that, in terms of 
classical mechanics, this entire C 2v symmetric subspace 
is forbidden for electrons. 

With perfectly symmetric escape excluded, we have to 
look for escape with almost symmetric electrons. Relax- 
ing the condition on the angular degrees of freedom, i.e. 
by allowing for arbitrary <fi — 4>i ~ 02 but p^ = 2p ( / )1 = 
—2pj, 2 (i.e. M = 0), we are led to consider the invariant 
subspace of C v symmetry (the C v subspace contains ob- 
viously the Civ subspace) Then the Hamiltonian of 
the system reads 



rr PR+P 2 Z , V* 



1 



VR 2 + Z 2 



2R 



(I) 



-2ZF(t). 



(4) 

Comparing (0| with we see that in the present 
case electrons can perform a symmetric bending motion 
around the field axis. In the C v subspace electrons may 
pass over the saddle considered previously and ionize but 
their momenta need not lie in the plane \<f)i — <\>i\ = it. 
In the quantum description there is no selection rule for 
electrons with fci = ki and with momenta that have 
the same projections on the polarization axis, but with 

101 - 02 | ^ 7T- 

The consequences of the selection rule should show up 
most clearly in angle resolved double ionization cross sec- 
tions, as they result in a zero for the symmetric momenta 
and an azimuthal angle difference of tt. However, it will 
also show up as a reduction in intensity in the angle in- 
tegrated cross sections, like the ones where only the par- 
allel momenta components of ionizing electrons are mea- 
sured P, Q . We will demonstrate this with calculations 
in 1+1-dimensional model below, but we emphasize that 
the cleanest verification is possible with fully momentum 
resolved cross sections as obtained in recent experiments 
1- 

Full three-dimensional (3D) calculations of the dou- 
ble ionization process are barely feasible [l^, H(| , so we 
take recourse to low-dimensional models. Aligned elec- 
tron models [lij are not suitable for our purpose since 
there the symmetric simultaneous escape of electrons is 
forbidden by the overestimated electron repulsion. The 
1+1-dimensional (1+1D) model proposed in [21J does al- 
low for symmetric escape by confining the electrons to 
move along one-dimensional (ID) tracks that form an- 
gles 7r/6 and 57r/6 with respect to the field polarization 
axis. The Hamiltonian of the system, in atomic units, 
reads 
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R2), 



(•5) 



where R\ and R 2 denote positions of electrons on the cho- 
sen tracks and pn t and pr 2 are conjugate momenta. In 
this model the rescattering process, single and sequential 
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double ionization are present but importantly also the 
symmetric simultaneous escape of electrons can be prop- 
erly simulated (HJ. Indeed, in the 3D case, there is a sin- 
gle saddle around which two electrons can escape simul- 
taneously and that saddle possesses two unstable direc- 
tions. The 1+1D model also has a saddle located at the 
same place as in the 3D case, R\ = R2 = — sgn[F(t)]R s , 
also with two unstable directions. Hence, the key topol- 
ogy elements of the phase space of the 3D case are prop- 
erly reproduced in the model. The symmetric subspace 
of the classical phase space in the model is defined by 
i?i = i?2 and pn t = Pr 2 . If, in quantum simulations, we 
start with the ground state of the unperturbed system 
symmetric electron escape is not forbidden. However, if 
we chose as an initial state the first excited state, which 
is antisymmetric with respect to electron exchange (in 
analogy to the 3 S state) the wavefunction vanishes at 
Ri = i?2- The symmetric subspace is then forbidden. 

We have performed numerical integration of the 
Schrodinger equation corresponding to the Hamiltonian 
([5]) with Coulomb singularities smoothed by the replace- 
ment 1/x — > + e with e = 0.6, starting with 
the ground and first excited states of the unperturbed 
Hamiltonian. The ground state energy, for the cho- 
sen parameter e is —2.83 while the first excited state 
is at —2.21. We have used a single cycle pulse of the 
form F(t) = Fo f (t) s'm(ujt) , where the pulse envelope 
f(t) = sin 2 (ujt/2), with u = 0.06 and F = 0.15. Sin- 
gle ionization in the case of the initial first excited state 
is very efficient (with probability equal to 0.98) and much 
more probable than in the case of the ground state (with 
probability 10~ 4 ). That is due to the higher energy of 
the first excited state, which allows electrons to easily 
pass over the Stark barrier and to escape near the peak 
intensity of the field. Probabilities for double ionization 
are comparable in both cases, about 10 -6 , thus indicat- 
ing that the difference in initial energy does not have an 
exessively strong effect on double ionization. 

Figure [1] shows electron momentum distributions in 
the double ionization of He for the initial ground and 
first excited states. The distributions have been ob- 
tained adopting the method proposed in [23| as described 
in [22^ . The momentum distributions are obtained by 
Fourier transforms of parts of the final wavefunctions in 
the regions |i?i|, |i?2| > 50 and they are Gaussian (of 
width 0.07) averaged in order to mimic the experimental 
resolution. For double ionization starting from the sym- 
metric ground state there is a symmetric distribution of 
parallel final momenta with a maximum along equal mo- 
menta, pRx = Pr 2 ■ That the distribution reaches only 
into the negative momentum sector is due to the short 
pulse which contains only one cycle, [24j |. When the ini- 
tial state is the antisymmetric state, the corresponding 
distribution shown in Fig. [TJd has a clear minimum along 
the axis pr 1 — pr 2 , reflecting the absence of the perfectly 
correlated escape, as predicted by the selection rule. The 
small residual contribution at the center is due to aver- 
aging with the Gaussian that models experimental reso- 
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FIG. 1: Electron momentum distributions for double ioniza- 
tion in strong linearly polarized laser field starting from the 
ground state of the unperturbed atom (a) and from first ex- 
cited state (b). The distributions are obtained from Fourier 
transforms of the parts of the final wavefunctions in the re- 
gions \Ri\, \Ra\ > 50 and are averaged with a Gaussian of 
width 0.07 in order to model experimental resolution. The 
time-dependance of the field is F(t) = Fof(t) sin(urf) where 
Fo = 0.15, uj — 0.06 and one-cycle pulse duration with an en- 
velope f(t) = sin 2 (tirf/2). The concentration of the outgoing 
momenta in one quadrant is a result of the single cycle laser 
pulse. The low intensity along the diagonal in panel (b) is the 
signature of the suppression of symmetric double escape by 
the selection. 

lution. 

The signatures of the effect can be enhanced, as shown 
in Fig. [2k., by considering the distributions of the mo- 
menta differences, pr 1 — pr 2 - One finds a strong max- 
imum at pr 1 — pr 2 = for double ionization starting 
from the symmetric ground state and a profound minium 
when the initial state is the antisymmetric first excited 
state. Note again that the minimum in Fig. [2Ji does not 
reach zero because of the averaging. But Fig. [5k clearly 
demonstrates that the resolution attainable in current ex- 
periments is sufficient to observe effects of the selection 
rule we analyze here. 

In Fig. [5Jd the distributions of momenta sum (equiva- 
lent to ion momentum distributions), i.e. pR x +Pr 2 , in 
double ionization of He atoms are plotted. We see that 
the momenta sum distribution shows no simple signature 
of the absence of highly correlated electron escape in the 
case of the initial first excited state except, possibly a 
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FIG. 2: The probability distributions for differences, pn 1 — 
Pr 2 , [panel (a)], and for sums, pr ± +Pr 2 , [panel (b)], of elec- 
tron momenta for double ionization in a strong linearly po- 
larized laser field. Data for initial ground state (red dashed 
curves) are compared with those obtained for the initial first 
excited state (solid black lines). The minimum in the differ- 
ence distribution in panel (a) reflects the suppresion of sym- 
metric escapes. All parameters are the same as in Fig. [TJ 

small narrowing of the distribution. By comparison, the 
existence of correlated pairs of electrons for the initial 
ground state case broadens the corresponding sum dis- 
tribution (dashed line). We also would like to point out 
that the asymmetry of the ion momentum distribution 
seen in Fig . reflects the carrier-envelope phase as ex- 
pected [241 ] . 



Finally, we would like to stress that the numerical re- 
sults presented are obtained in the reduced dimension- 
ality model where the ID electron tracks are chosen in 
two-dimensional space containing the polarization axis. 
Direct application of the reduced 1+1D model to exper- 
imental data is thus possible only if events with electron 
momenta lying in a plane that contains the field polar- 
ization axis are collected only. 

We have shown that the selection rules of Maulbetsch 
and Briggs [11] imply a difference in the final momenta 
distributions in double ionizatoin by strong linearly po- 
larized laser pulses depending on the symmetries of the 
initial state. For double ionization in He starting from 
the metastable state 3 S correlated electron escape with 
final electron momenta symmetric with respect to the 
polarization axis is suppressed. Numerical simulations 
in a reduced dimensionality model show a clear dip in 
the momentum distributions, and demonstrate that the 
effect appears despite differences in the initial energy. 
The effect should be within reach of todays experiments 
with their resolution and detailed kinematic character- 
ization of the final state. Extensions to double ioniza- 
tion in molecules are immediate. A prime candidate are 
O2 molecules with triplet initial state oriented along the 
polarization axis. Recent experiments with aligned N2 
molecules JIt} probed a singlet initial state only. 
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